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application  in  Imaging  systems.  Computer  simulation  studies  were  made 
for  coherent,  partially  coherent,  and  incoherent  imaging.  These  were 
done  using  various  conditions,  including  variable  signal-to-noise  ratio 
and  phase  aberration,  but  with  monochromatic  illumination.  We  showed 
-that,  for  incoherent  imaging,  a large  degree  of  thinning  down  to  6% 
of  the  number  of  full  array  elements  (4096)  produces  image  quality 
comparable  to  that  produced  by  the  full  array.  For  broad-band 
polychromatic  incoherent  imaging, even  larger  degrees  of  thinning, down 
to  3Z,  were  similarly  effective.  The  fractional  degree  of  thinning 
allowable  for  incoherent  imaging  in  principle  is  inversely  proportional 
to  the  square  root  of  the  number  of  elements  in  the  full  array,  which 
number  is  equivalent  to  the  number  of  elements  resolved  in  the  full 
array.  VThe  arrays,  which  are  easy  to  design,  degrade  gracefully 
when  large fractions  of  the  elements  are  removed.  For  coherent  imaging 
application^.,  random  arrays  do  not  offer  significant  advantages  over 
full  arrays.  N. 
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SECTION  1 


INTRODUCTION  AND  SUMMARY 

Recent  requirements  for  very  high  resolution  imaging  systems  have 
produced  a great  demand  for  apertures  with  high  resolution,  low  sidelobe 
levels,  and  high  gain.  To  achieve  the  resolutions  desired,  extremely 
large  apertures  are  needed.  The  realization  of  such  apertures  is  gen- 
erally considered  to  be  either  economically  or  technically  infeasible. 
Table  1 gives  some  representative  numbers  for  such  arrays.  The  .lumbers 
are  astronomically  large  in  many  cases.  Hence,  attention  has  been 
focused  on  large  arrays.  Arrays  with  uniformly  spaced  elements  produce 
undesirable  grating  lobes.  For  these  arrays,  there  exists  a maximum 
spacing  between  elements  if  the  grating  lobes  are  to  be  removed  from 
the  object.  The  number  of  elements  required  in  such  uniform  arrays  would 
be  very  large  and  in  direct  proportion  to  the  aperture  dimension. 
Therefore,  the  investigation  of  arrays  with  a reduced  number  of  ele- 
ments (the  "thinned  arrays"),  which  possess  the  desired  characteristics 
of  the  full  arrays,  has  been  of  recent  Interest, 

Arrays  with  incommensurable  element  spacings  usually  have  aperiodic 
radiation  pattern  functions.  Consequently,  the  number  of  elements 
required  is  not  directly  determined  by  the  aperture  dimension.  However, 
there  is  no  general  theory  available  for  the  algorithmic  design  of  this 
class  of  arrays.  Many  of  the  designs  to  date  have  been  based  on  trial- 
and-error  methods  using  high-speed  computers.  But  there  is  never  any 
assurance  that  any  particular  trial  will  produce  a successful  design. 

One  class  of  arrays  exists  that  possesses  the  desired  character- 
istics of  high-resolution  imaging  systems  and  yet  requires  surprisingly 
few  elements  to  achieve  these  features.  This  class,  random  arrays,  is 
well  established  in  antenna  theory  and  has  been  generally  applied  to 
signal-detection  applications.  The  properties  of  such  arrays  are 
approached  from  a probabilistic  point  of  view,  even  though  the  per- 
formance of  any  particular  array  itself  is  not  probabilistic.  That  is, 
once  the  positions  of  the  array  elements  are  determined,  the  properties 
of  the  random  array  are  completely  deterministic.  A prime  advantage 


of  these  arrays  13  that  predictions  on  the  array's  performance  can  be 
made  before  carrying  out  detailed  computations.  Therefore,  designing 
the  array  is  reduced  to  playing  a game  in  which  the  probability  of 
success  is  determined  to  be  sufficiently  large  before  any  actual  evalu- 
ation of  the  final  array  design  is  attempted.  Other  advantages  of 
random  arrays  are  discussed  below. 

For  random  arrays,  the  sidelobe  level  is  closely  related  to  the 
number  of  elements  and,  to  a much  lesser  degree,  to  the  total  aperture 
dimension.  The  resolution,  or  main  beam  width,  depends  on  the  total 
aperture  dimension;  the  directive  gain  is  proportional  to  the  number  of 
elements  used  if  the  average  spacing  is  large.  Consequently,  extremely 
high  resolution  can  be  achieved  with  very  few  elements.  On  the  other 
hand,  for  a given  number  of  elements,  higher  resolution  can  be  obtained 
by  spreading  these  elements  over  a large  aperture,  with  the  sidelobe 
levels  remaining  substantually  the  same  and  the  directive  gain  constant. 

We  have  applied  the  concept  of  random  arrays  to  imaging  systems  to 
determine  their  capabilities  and  limitations  for  imaging  applications. 

To  the  best  of  our  knowledge,  this  is  the  first  demonstration  of  the 
utilization  of  thinned  random  arrays  in  Imaging  systems.  Computer 
simulation  studies  have  been  carried  out  here  for  incoherent,  partially 
coherent,  and  coherent  imaging.  These  investigations  were  performed 
using  various  conditions , including  variable  signal-to-noise  ratio  (S/N) 
and  phase  aberration  in  the  Imaging  system,  with  monochromatic  and  poly- 
chromatic illumination  of  the  object.  However,  no  attempt  was  made  to 
post-process  the  reconstructed  images  for  any  type  of  enhancement. 

The  results  show  that  good  quality  images  can  be  obtained  for 
ranuom  arrays  with  very  few  elements,  even  down  to  only  6%  of  the 
number  of  elements  of  a uniform  full  array  of  4096  elements.  (The 
number  of  elements  in  the  full  array  is  equivalent  to  the  number  of 
elements  resolved  by  the  full  array  in  its  field  of  view  (FOV)).  These 
results  hold  for  a low  (e.g.,  > 3)  S/N  and  a high-intensity  background 
in  the  image.  For  broadband  polychromatic  incoherent  imaging,  the 
number  of  elements  can  be  even  further  reduced.  In  addition,  this 
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Image  quality  can  be  achieved  by  many  distinct  random-array  configurations. 
The  image  resolution  can  be  maintained  as  the  number  of  array  elements 
is  reduced  if  the  aperture  dimension  is  fixed  for  each  random  thinned 
array.  The  quality  of  the  image  degrades,  as  shown  by  the  appearance 
of  higher  (random)  background  levels  as  the  number  of  array  elements 
is  reduced.  This  is  a manifestation  of  the  increased  sldelobe  level  of 
the  random  array  as  the  number  of  elements  is  reduced.  For  a given 
number  of  elements,  the  simulation  results  exhibited  an  improvement 
in  image  resolution  as  the  number  of  elements  was  spread  over  larger 
aperture  dimensions.  The  overall  background  levels  in  these  images 
remained  essentially  constant.  These  encouraging  results,  whlbh  were 
obtained  without  utilizing  any  averaging  or  interpolation  post-processing 
techniques  in  the  aperture  plane  for  image  enhancement,  lead  us  to 
believe  that  similarly  Impressive  performance  can  be  achieved  in 
practice. 

We  have  studied  the  effects  of  background,  S/N,  and  phase  aberra- 
tions (such  as  those  caused  by  turbulence  and  any  figure  disposition 
in  the  case  of  incoherent  imaging).  For  coherent  imaging,  the  effects 
of  speckle  were  examined,  and  a method  of  speckle  elimination  was 
studied.  The  determination  of  the  dependence  of  the  maximum  thinning 
factor  possible  to  retain  good  quality  in  the  image  was  investigated 
as  a function  of  the  pertinent  parameters,  Including  the  totel  number 
of  elements  and  S/N. 

To  present  a comprehensive  exposition  of  our  results,  we  have 
incorporated  in  this  report  the  principle  results  of  the  previous 
reporting  periods  along  with  the  two  studies  performed  in  this  period. 

In  this  year,  a study  was  made  of  the  use  of  broadband  polychromatic 
radiation  incoherent  imaging.  Partially  coherent  imaging  was  modeled 
as  well,  and  computer  simulations  were  performed. 
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SECTION  2 


TECHNICAL  DISCUSSION 


A.  IMAGING  CONCEPTS 

The  imaging  process  may  be  viewed  in  two  complementary  ways. 

Although  mathematically  equivalent  and  related  by  theorems  of  Fourier 
transformation,  one  view  or  the  other  reveals  special  insights  and 
simplifications,  depending  on  the  particular  circumstances.  We  may, 
for  example,  view  the  imaging  process  as  a convolution  of  the  object 
field  and  the  radiation  or  antenna  pattern  of  the  aperture,  as  is  the 
case  for  coherent  imaging.  If  the  aperture  is  a large  lens  or  collect- 
ing mirror,  then  its  radiation  or  diffraction  pattern,  which  is  the 
Fourier  transform  of  the  aperture  configuration,  will  be  a sharply 
peaked  function  in  angular  extent  with  only  small  subsidiary  sidelobes. 
Convolving  this  pattern  with  the  object  field  produces  an  image  field 
the  quality  of  which  depends  on  the  strength  and  degree  of  sharpness 
of  the  mainxobe  and  the  weakness  and  disposition  of  the  sidelobes. 

Alternatively,  we  may  view  the  imaging  process  by  considering  the 
Fourier  transformation  of  the  object  into  its  spectrum  of  spatial  fre- 
quencies. If  only  the  intensities  of  the. image  and  object  are  of 
interest,  then  the  Fourier  transformation  of  the  intensity  representa- 
tion of  the  object  will  appear  in  the  far  field  at  the  aperture,  or 
entrance  pupil,  of  the  imaging  instrument  as  the  (complex)  visibility 
pattern;  it  is  sometimes  called  the  partial  coherence  transform  vis-a-vis 
the  Van  Citert-Zernike  theorem.*  This  coherence  transform  is  sampled 
to  some  degree  in  the  aperture  plane,  depending  on  the  nature  and  size 
of  the  collecting  aperture,  Th  t amplitude  of  the  transform  is  the 
classical  visibility  of  fringes  as  a function  of  spatial  frequency  given 
by  the  separation  of  the  aperture  elements.  The  phase  of  the  trans- 
form is  the  spatiul  phase  of  these  fringes  referred  to  the  optical  axis 
or  focus  of  the  system.  This  sampled  function  is  then  inversely 

*This  theorem  is  the  spatial  analog  of  the  Wiener-Khintchine  theorem 
relating  the  power  spectrum  and  temporal  correlation  for  stationary 
processes. 

13 


transformed  to  form  the  image.  The  quality  of  the  image,  from  this 
point  of  view,  depends  on  the  degree  to  which  the  irradiance  of  the 
object  has  been  sampled  in  the  aperture  plane  and  the  disposition 
of  the  sampling  elements.  Insufficient  and  regular  sampling  of  the 
aperture  will  produce  an  image  not  only  lacking  in  all  the  correct 
details  at  the  missing  (i.e.,  unsampled)  spatial  frequencies,  but  will 
also  produce  vivid  spurious  image  information  caused  by  large  uncom- 
pensated Fourier  terms. 

To  amplify  the  description  of  imaging  and  to  cast  these  concepts 
in  more  concrete  terms,  it  is  instructive  to  put  the  results  in  more 
formal  mathematical  terms.  The  mathematical  representations  of  the 
imaging  process  take  especially  simple  form  in  the  limiting  cases  of 
completely  coherent  and  completely  incoherent  imaging.  These  cases 
are  of  practical  interest  as  well,  for  they  are  good  approximations 
to  conditions  encountered  most  often.  However,  we  did  not  limit  our 
study  to  these  cases. 


1.  Coherent  Imaging 

Coherent  imaging  systems  are  linear  in  complex  field  amplitude. 
From  diffraction  considerations,  the  field  distributions  of  the  image, 
U^(x,y) , and  object,  Uo(u,v),  are  related  by  the  convolution  Integral 


l^Cx.y)  - JJ  h(x-u,  y-v)  Uq(u,v)  dudv  , (1) 

obj 

where  h(x-u,  y-v)  is  the  impulse  response  of  the  imaging  system.  A 
generalized  schematic  of  an  imaging  system  is  shown  in  Figure  1.  The 
pupil  function  P(x',y’)  has  the  property  that 


P(x\y’) 


1 within  the  aperture  configuration 
0 elsewhere. 


(2) 


This  function  accounts  for  the  finite  extent  of  the  pupil,  or  aperture, 
of  the  imaging  system.  The  system  impulse  function  is  directly  related 
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to  the  pupil  function  by  a Fourier  transformation: 


h(x,y)  - //  P(a,e)e“j2lT<xa+ye)  dctdb  . (3) 

aper 

Thus,  for  a given  wavelength,  the  Imaging  properties  of  a system  are  com- 
pletely determined  by  the  pupil  function,  or  aperture,  of  the  imaging 
system. 

Two  quadratic  phase  terms  that  normally  appear  in  Eq.  3 have  been 
neglected  in  the  diffraction  formulation  of  the  imaging  system.  These 
terms  convey  phase  curvature  information  and  can  be  dropped  directly  if 
imaging  is  considered  to  *ie  between  two  spherical  surfaces.  For  plane 
surfaces,  these  terms  can  also  be  neglected.  The  justification  is 
based  on  the  fact  that  the  resultant  image  is  usually  sensed  with 
detectors  that  respond  to  intensities.  One  of  these  quadratic  phase 
terms  only  modifies  the  phase  distribution  of  light  and  cannot  affect 
the  resultant  detected  image  intensities.  The  second  quadratic  phase 
term  can  be  neglected  also  if  the  imaging  system  maps  neighboring  points 
in  the  object  plane  onto  neighboring  points  in  the  image  plane.  This 
condition  is  satisfied  for  most  cases  where  the  object  is  of  small 
angular  extent.  Mathematical  details  of  these  concepts  can  be  found  in 
J.W.  Goodman,  Introduction  to  Fourier  Optics. 

400C - t 


Figure  1.  Schematic  of  generalized  imaging  system. 
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The  convolution  integral,  Eq.  1,  can  be  conveniently  recast  into 
its  spatial  frequency  spectrum  representation  by  Fourier  transformation. 
Denoting  the  Fourier  transforms  of  U and  h by  0 and  K,  respectively, 
the  Fourier  representation  of  Eq.  1 is 


u.(f  »f  ) 

i x y 


h(f  ,f  > U (f  ,f  ) 
x y o x y 


(4) 


where  f and  f are  the  Fourier  spatial  frequency  components.  Here, 

^ y 

h(f  ,f  ) is  known  as  the  coherent  transfer  function  (CTF).  Recalling 
x y 

from  Eq.  3 that  h(x,y)  is  already  the  Fourier  transform  of  the  pupil 
function  P(x,y),  then 


h<vy 


p<vy 


(5) 


and  Eq.  4 becomes 


wy  ■ p<vy  vf*<y  • 


(6) 


where  magnification  and  image  inversion  have  been  ignored.  The  rela- 
tionship  between  the  coherent  transfer  function  h and  the  pupil  func- 
tion (or  aperture)  P is  a direct  and  simple  one.  In  the  Fourier 
transform  representation,  therefore,  the  pupil  function  (or  aperture) 
samples  the  Fourier  spatial  frequency  components  of  the  object  trans- 
form field.  The  image  can  be  reconstructed  from  the  sampled  frequency 

components  by  an  inverse  transformation  of  U.(f  ,f  ). 

l x y 

As  Eq.  6 shows,  diffraction  effects  of  the  imaging  system  are 
manifested  by  the  finite  extent  of  the  pupil  function.  That  is,  high 
spatial  frequency  components  of  the  object  transform  field  are  "cut 
off"  or  not  sampled.  From  another  point  of  view,  the  high-frequency 
components  of  the  object  field  diffract  to  a greater  angular  extent 
than  the  low-frequency  components:  they  are  not  "intercepted"  by  the 
finite  aperture  of  the  imaging  system.  To  reconstruct  only  the  high- 
frequency  content  of  the  object  (e.g. , the  edge  content),  it  would  be 
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necessary  to  sample  only  the  high  spatial  frequency  components  of  the 
object  transform  field.  In  this  case,  an  annular  aperture  would  be 
adequate. 


2.  Incoherent  Inaging 

Incoherent  imaging  systems  can  be  viewed  from  the  same  approach  as 
described  above.  For  such  incoherently  illuminated  objects,  the  imag- 
ing system  is  linear  in  intensity  distribution.  From  diffraction 
considerations,  and  representing  the  image  and  object  field  intensities 
by  I^(x,y)  and  Iq(u,v),  respectively,  these  field  intensities  are 
related  by  the  convolution  integral 


|h(x-u,  y-v) j ‘ 


IQ(u,v)  dudv 


(7) 


obj 

where  h(x-u,  y-v)  is  again  the  impulse  response  of  the  imaging  system. 
By  forming  the  intensity  of  Eq.  1 and  imposing  the  condition  of  statis- 
tical independence  (incoherence)  on  the  phasor  amplitudes  across  the 
object,  Eq;  7 can  be  derived  from  Eq.  1.  In  general,  the  statistical 
nature  (or  coherence)  of  the  object  phasor  amplitude  determines  the 
form  of  the  imaging  equations. 

The  image  field  intensity  is  also  determined  by  the  system  impulse 

response  h(x,y),  which,  in  turn,  is  related  to  the  pupil  function 

P(x',y')  by  Eq.  3.  However,  the  mapping  in  this  case  is  carried  out 
2 

by  |h(a,8) | , the  magnitude  squared  of  the  system  impulse  function. 

2 

The  function  |h|  , known  as  the  point  spread  function  (PSF),  gives, 
for  finite  apertures,  the  degradation  of  a 6-function  (point)  source 
through  the  imaging  system. 

Because  of  the  convolution  fora  of  Eq.  7,  it  can  be  recast  into 
its  spatial  frequency  spectrum  representation  aleo.  Denoting  the 

i 1 2 *•»* 

Fourier  transforms  of  I and  |h|  by  I and  H,  respectively,  the  Fourier 
representation  of  Eq.  7 is 


I,(f  ,f  ) = 8(f  ,f  ) I (f  ,f  ) 
i x y x’  y o x y 


(3) 
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The  function  H(f  , f ) is  commonly  called  the  optical  transfer  function 
x y 

(OTF)  of  the  incoherent  imaging  system,  and  its  modulus  |H(f  ,f  )|  is 

x y 

known  as  the  modulation  transfer  function  (MTF).  In  the  Fourier  trans- 
form representation,  the  OTF  acts  as  a complex  weighting  factor  in 
mapping  the  Fourier  spatial  frequency  components  of  the  object  field 
intensity  into  image  Fourier  spatial  frequency  components.  The  image 
can  thus  be  reconstructed  from  the  mapped  frequency  components  in 

transform  space  by  an  inverse  transformation  of  I.  (f  ,f  ) 

i x y 

Both  the  CTF,  h(f  ,f  ),  and  OTF,  H(f  ,f  ),  involve  the  system 
x y : x y 

impulse  response  h(a,3).  There  exists  a specific  relationship  between 
them.  The  autocorrelation  theorem  can  be  used  to  show  that 


H(f  ,f  ) - 
x’  y 


//  Kfa, 


3)  h (a+f  , 3+f  ) dad 3 
x y 


and  thus  that  the  OTF  is  the  autocorrelation  of  the  CTF,  Using  Eq.  5, 
the  above  relationship  can  be  reduced  further  to  show  that  the  OTF 
H(fx,fy)  is  the  autocorrelation  of  the  pupil  function  (or  aperture) 
PCx’.y'): 


B(fx,f  ) = //  r'<*,3)  P*(a+fx,  3+f) 


There  is  an  interesting  and  important  geometric  interpretation  to  the 
autocorrelation  function  of  the  aperture  distribution.  It  counts  or 
measures  the  distribution  of  all  element  pairs,  discrete  or  differen- 
tial, as  a function  of  pair  separations,  independent  of  the  absolute 
position  of  the  elements.  This  may  be  readily  seen  by  considering  a 
one-dimensional  array  of  discrete  elements  irregularly  disposed  on  a 
regular  lattice.  The  first  point  (zero  separation)  of  the  autocorrela- 
tion function  is  the  result  of  multiplying  the  array  with  itself,  and 
merely  counts  the  number  of  elements  in  the  array.  Displacing  the 
array  over  itself  by  one  unit  produces  overlaps  wherever  there  were 
unit  separations  in  the  array  itself;  on  multiplication,  the  number  of 
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pairs  of  unit  separation  is  counted.  Similarly,  the  larger  displacements 
measure  the  existence  of  larger  separation  pairs.  Since  the  separation 
of  a pair  exists  in  spatial  frequency  space,  the  OTF  for  incoherent  light 
can  be  viewed  as  a weighting  function  with  a modulus  equal  to  the  number 
of  pairs  of  elements  at  each  spatial  frequency. 


Fartiallv  Coherent  Imaging  Transfer  Function 


The  concept  of  a linear  system  impulse  response  can  still  be  used, 
even  in  the  case  of  partially  coherent  imaging,  although  its  representa- 
tion is  somewhat  more  complex.  Recalling  the  definition  of  the  mutual 
coherence  function  r(x^,X2»t^t2)  as  the  cross  correlation  of  the  fields 
at  U(x^,t^)  and  U(x2,t2),  and  assuming  the  fields  are  stationary  in 
time,  T may  be  written  as  a function  of  difference  in  time  T * t^  - t2: 


r(x1,x2,T)  » <U(xr  t-K)  U (x2,t)>  , 


where  the  bar  notation  (xi>  indicates  a two-dimensional  spatial  vector 
(x  ,y  ).  The  coherence  function  satisfies  the  wave  equation  and 
therefore  propagates,  as  does  light  itself,  with  the  effects  of  diff ration. 
As  a result,  the  Fourier  components  of  the  coherence  function  for  each 
carrier  frequency  at  the  image  may  be  related  to  that  at  the  object  by 
an  impulse  response  function  not  unlike  that  for  fields  themselves: 


’ f ’ ’ * — __  

r(x1,x2,a))  = h (Xj-u^,  w)  h (x2-u2,  w)  r(uj,u2,tt))  du-jd^  , (12) 


obj ect 


where  r(w)  is  the  Fourier  transform  of  T(t).  I ] 

Assuming  the  quasi-monochromatic  approximation  j j 

t 1 

) E 
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r(u1,u2,a))  - f(u1u2,0)  S(w  - 0))  , (13)  I | 


we  may  write,  after  retransformation  and  simplification  of  the  integral. 


r(xr 

where  to  is  the  mean  carrier  frequency.  This  relates  the  mutual  coherence 
of  the  object  and  image.  Since  the  intensity  in  the  image  equals 
T(x1,x1,0)  of  the  image,  we  may  write 

KXf)  » j jj  [h(x1~u1,  to)  h*(x2-u2,  to)]  r(U;L,u2,0)  duxdu2  , (15) 

obj ect 

where  the  bracketed  quantity  may  be  considered  to  be  the  impulse  response. 
For  numerical  or  experimental  analysis,  the  image  intensity  can  be  generated 
from  Eq.  15.  Equivalently,  for  a given  form  of  T in  terms  of  the  field 
or  intensity,  Eq.  15  can  be  written  entirely  as  a linear  superposition 
of  object  mutual  intensities.  r(u^,u2»0)  is  sometimes  called  the  mutual 
intensity  J(u^,u2>. 

Particular  explicit  models  of  the  coherence  F(u^,u2>  of  the  object 
fields  allow  particular  development  of  the  image  intensity  I(x^)  in 
terms  of  the  object  fields  or  intensities.  Of  course,  Fourier  trans- 
forming the  impulse  response  relation  (Eq.  15)  by  the  convolution 
theorem  recasts  it  into  the  form  of  a linear  transfer  function.  Four- 
dimensional  convolution  or  Fourier  transforms  on  sizable  images  will 
exceed  the  capabilities  of  most  computers. 

To  simplify  our  computer  calculations,  we  consider  an  optical  system 
with  transmission  function  for  mutual  intensity  K(u^,x^).  Using  scale- 
normalized  coordinates  where  the  object  point  and  its  Gaussian  image 
have  the  same  coordinate  numbers,  the  transmission  of  the  mutual  intensity 
is  given  by: 
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(16) 
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A second  major  assumption  we  will  make,  besides  monochromaticity,  is 
that  the  object  is  small  enough  to  form  an  isoplanatic  region.  In  this 
case,  the  transmission  function  K is  a function  only  of  the  differences 
in  coordinates,  for  example,  K(x  - u),  and  not  of  the  coordinates  them- 
selves. Then  Eq.  16  becomes 


Jr  (x1,x2)  - | ||  J0(ultu2)  K(x1-u1)  K*(x2-u2)  dujdujj  . (17) 

object 

This  relationship  is  a convolution  and  we  can  apply  the  convolution 
theorem  to  the  Fourier  transforms  of  the  parameters  to  obtain 

J<f..g;f,.g'>  - J(f,g;f’,g')  K(f,g)  K*(-£,,-g,>  . (18) 

The  transmission  of  mutual  intensity  is  seen  to  still  be  a four- 
dimensional  linear  system.  Note  also  that  the  Fourier  transforms  are 
four  dimensional. 

The  *>.bove  formulation  is  simple  in  form,  but  still  leads  to  signifi- 
cant problems  in  terms  of  implementation  in  digital  image  processing. 

For  example,  in  our  simulations  the  smallest  basic  image  is  represented  by 

an  array  of  128  x 128  pixels.  To  represent  the  mutual  intensity  for  this 
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image  would  require  2 complex  numbers,  which  is  far  in  excess  of 
available  storage.  Also,  computing  the  four-dimensional  Fourier  transforms 
for  such  an  image  is  computationally  prohibitive,  even  using  fast  Fourier 
transform  techniques. 

To  mitigate  the  data  storage  problem,  we  use  the  degree  of 
coherence,  y,  defined  by 


_ __  •Ku1»u2) 

“<ui'u2)  - TmpTmj 


(19) 


To  understand  how  this  parameter  can  help  to  simplify  the  problem, 
consider  the  case  of  partial  coherence  aiising  from  a uniform  extended 
incoherent  illumination  source.  By  the  VanCittert-Zernike  theorem,  if 
the  distance  between  the  points  (u^)  and  (u2)  is  small  compared  to  the 
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distance  to  the  source,  the  degree  of  coherence  la  given  by  the  normalized 
Fourier  transform  of  the  source  intensity  distribution.  Thus,  for  a 
Gaussian  source,  y will  be  Gaussian.  For  a square  source,  it  will  be  the 
two-dimensional  sine  functions;  for  a circular  source,  it  will  be  the 
Bessel  function  of  first  kind  and  first  order. 

Basically  then,  the  degree  of  coherence  of  the  illumination  on  the 
object  will  be  a function  only  of  the  difference  in  coordinates  for  those 
cases  of  interest.  Thus  the  mutual  intensity  of  the  object  can  be 
written 


VvV  " VU1  - u2)  /f0<u2T  ' (20) 

That  is,  we  can  represent  the  mutual  intensi*  lea  by  two  two-dimensional 
functions:  y and  I. 

The  problem  of  computing  four-dimensional  FFTs  remains.  If  Eq.  20 
is  examined,  we  see  that  this*  representation  has  terms  that  are  separable 
by  points  (the  intensity  function),  and  it  has  a term  y that  is  normally 
separable  by  coordinates  (as,  for  example,  with  Gaussian  or  sine  forms  for 
the  degree  of  coherence  of  the  illumination) . If  the  degree  of  coherence 
would  be  approximated  by  a function  that  were  totally  separable  In  all 
coordinates,  then  Eq. . 20  could  be  separated  by  points  into  a product  of 
two  two-dimensional  Fourier  transforms  which  would  then  split  into  a product 
of  two  two-dimensional  transforms,  which  can  be  computed  easily.  Extending 
this  approach,  we  looked  for  an  approximation  for  the  modulus  of  y that  was 
a Bum  of  totally  separable  terms  and  was  a sum  only  of  the  absolute  values 
of  the  difference  in  terms.  To  this  we  added  a complex  phase  that  was 
separable  by  points,  as  are  the  intensity  terms  in  Eq.  20.  In  each  axis, 
the  functions  we  selected  for  a cosine  series  representation  of  y are  of 
the  form 

2 An  C0S  ^ (W  * 

n=0 


(21) 


where  a is  a parameter  that  determines  the  coherence  radius.  For  example, 
if  a is  zero,  the  functions  yield  coherent  imaging.  Many  realistic 
coherence  functions  can  be  well  represented  by  such  a cosine  representation. 

Note  that  this  function  is  separable  by  expanding  the  terms  using 
trigonometric  identities  and  that  it  is  symmetric  in  and  u2*  This 
results  in  a representation  for  the  mutual  coherence  of  the  form 

_ _ Z _ _ (u. ) i<t>tu,)  

Jo<Ul‘U2)  " Z W W 6 e * (22) 

n*l 

Since  the  optical  transmission  function  also  separates,  as  it  is  a set  of 
ones  and  zeros  distributed  over  the  aperture,  the  entire  transmission  as 
given  by  Eq.  18  also  separates  into  a sum  of  separable  functions. 

Appendix  C gives  an  example  of  a six-term  cosine  series  representa- 
tion of  y (Eq.  21)  and  of  the  resulting  121-term  separable  representation 
(Eq.  22).  The  result  is  plotted  in  Figure  2. 

B.  THINNED  ARRAY  REQUIREMENTS  FOR  IMAGING 

The  two  viewpoints  of  imaging  discussed  above  (i.e.,  the  convolu- 
tion integral  and  transfer  function  representations)  permit  a revealing 
approach  for  the  application  of  thinned  arrays  to  imaging  systems.  To 
achieve  high  resolution  in  the  linage,  the  high  spatial  frequency  content 
of  the  object  must  be  maintained.  Referring  to  Eq.  6 for  the  coherent 
case,  we  note  that  a large  aperture  is  required  to  sample  the  high 
spatial  frequency  content  of  the  object  transform  field.  Relating  this 
requirement  to  the  convolution  integral  representation  (Eq.  1)  and 
recalling  the  relationship  of  Eq.  3 between  the  system  impulse  function 
h(a,B)  and  the  aperture  P(x',y'),  it  is  seen  that  a large  aperture 
implies  a sharply  peaked  system  impulse  function  with  low  subsidiary 
sidelobe  levels.  *In  this  case,  the  contribution  to  the  image  field 
intensity  distribution  will  be  given  mainly  by  the  sharply  peaked  main- 
lobe  of  the  aperture,  with  relatively  small  contributions  coming  from 
the  low-level  sidelobes  in  the  convolution  process. 

For  the  incoherent  case,  similar  interpretations  can  be  given, 
although  the  aperture  (or  pupil  function)  P(x’,y')  appears  now  in  a less 
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direct  way  through  the  system  impulse  function  |h(a,3)|  , which  repre- 
sents the  received  power  distribution  of  the  array  pattern.  Again, 
referring  to  the  transfer  function  representation  (Eq.  8)  and  recalling 

the  autocorrelation  relation  of  the  aperture  to  the  OTF  H(f  ,f  ),  Eq.  10, 

x , y 

it  is  seen  that  high-resolution  images  require  detection  at  large  sepa- 
rations in  the  spatial  frequency  plane.  This  translates  again  into  a 
large  aperture.  Viewed  from  the  convolution  integral  diffraction  point 
of  view, Eq.  7,  the  system  impulse  function  (or  point  spread  function) 

J h | must  possess  a sharply  peaked  mainlobe  with  relatively  low  level 
sidelobes.  In  such  a case,  the  subsidiary  sidelobe  contributions  will 

not  be  sigi  leant  in  the  convolution  integral  and  therefore  will  not 

* 

appreciably  degrade  the  image. 

As  pointed  out  in  Section  1,  high-resolution  imaging  systems  may 
be  infeasible  because  of  the  large  apertures  required.  The  necessity 
of  utilizing  a much  smaller  number  of  elements  to  accomplish  the 
requirements  of  the  full  aperture  is  obvious.  Stated  most  simply,  the 
task  is  to  design  thinned  arrays  that  possess  the  desired  characteristics 
of  large  apertures  for  high-resolution  Imaging  applications:  sharply 
peaked  mainlobe,  low-level  sidelobes,  and  high  directive  gain. 

One  obvious  way  of  reducing  the  number  of  elements  required  by  a 
full  array  while  still  recovering  the  total  information  content  of  the 
object  is  to  impose  the  sampling  theorem.  Sampling  in  the  aperture 
plane  with  a periodic  array  of  receiver  elements  reproduces  an  image 
of  the  object  that  is  replicated  at  a distance  inversely  proportional 
to  the  sampling  spacing.  Spacing  the  sampling  locations  at  the  reci- 
procal of  the  object  dimension  and  scaling  appropriately  for  wavelength 
and  target  range  cause  the  replicas  to  just  touch.  This  is  the  sampling 
interval  dictated  by  the  sampling  theorem.  Spacing  at  finer  intervals 
spreads  the  replicas  farther  apart,  thereby  reducing  the  so-called 
aliasing  effect.  However,  this  is  achieved  at  the  expense  of  requiring 
more  array  elements,  while  not  increasing  the  information  content  of 
the  image.  By  truncating  the  sampled  space,  each  replica  has  the  edge 
ringing  that  even  a full  diffraction  limited  finite  aperture  would 
suffer. 
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These  properties  can  be  demonstrated  by  defining  the  convenient 
sampling  function  commonly  called  the  comb(x,a),  which  is  |a|E6(x  - na)  . 
The  comb  samples  a function  at  the  values  x ■ na  when  it  is  multiplied 
by  that  function.  In  addition,  convolving  the  comb  with  a function  f(x) 
replicates  that  function  at  the  interval  a: 


(1/a) comb  (x,a)  © f(x)  * Ef(x  - na) , 


where,  for  convenience,  © represents  convolution.  The  replicas  may 
overlap  depending  on  the  extent  of  the  function  compared  to  the  inter- 
val a.  The  Fourier  transform  of  a comb  is  another  (reciprocal)  comb. 

As  may  be  recalled  from  the  relationship  between  grating  periodicity 
and  grating  lobes  or  orders  in  diffraction  theory,  the  Fourier  transform 
of  comb  (x,a)  is 


comb 


I6(x'  - n/a) 


By  the  convolution  theorem,  if  a sampled  function  is  Fourier  transformed, 
then  a replicated  transform  of  the  function  is  produced: 


Fourier  transform  of 


comb 


f(x’) 


f(x)  © comb (x, a)  . 


Because  of  the  reciprocal  nature  of  the  transformed  comb,  the  closer 
the  original  samples,  the  farther  apart  are  the  replicas  in  the  convolu- 
ted transform.  If  an  object  f(x)  is  bounded  at  xnjax»  then  its  replicas 
will  just  touch  if  the  object  transform  f(x')  is  sampled  at  a * 1/x^^. 
Finer  sampling  of  the  object's  transform  f(x')  spreads  the  replicas 
apart.  This  may  be  viewed,  in  the  case  of  the  imaging  problem,  as 
spreading  apart  the  grating  lobes  (the  convolution  comb)  of  the  aperture 
array  (the  sampling  comb) . Limiting  the  required  infinite  number  of 
samples  needed  for  exact  reproduction  of  each  replica  to  a sample  of 
finite  extent  distorts  the  reproduction  in  exactly  the  same  way  as  does 
a truncated  full  aperture  (i.e.,  continuous  sampling).  It  produces  all 
the  well  known  distortions  anu  edge  "ringing"  of  diffraction  by  a trun- 
cated aperture.  The  number  of  sampling  points  required  equals  the 
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number  of  resolution  elements  seen  on  the  object.  For  the  coherent  case, 
the  number  of  sampling  points  is,  in  fact,  the  number  of  elements  in  the 
array;  for  the  incoherent  case,  each  sampling  point  is  derived  from  a 
pair  of  elements.  N elements  provide  N(N-l)/?  pairs  including  possible 
redundancies. 

In  determining  the  best  thinned  array  configurations  for  high- 
resolution  imaging,  the  vast  amount  of  established  analysis  of  antenna 
array  theory  can  be  fruitfully  applied  to  imaging  systems.  In  particu- 
lar, the  concept  of  random  arrays  appears  to  be  an  optimum  approach  to 
this  problem.  Although  random  arrays  have  been  analyzed  in  detail  for 
use  in  microwave  antenna  theory,  their  utilization  has  been  restricted 
mainly  to  signal  detection  in  the  microwave  region.  To  the  best  of  our 
knowledge,  their  application  to  imaging  systems  has  not  been  explored. 

The  properties  of  random  arrays  and  their  superiority  for  high- 
resolution  imaging  applications  are  discussed  in  the  next  section. 

C.  THINNED  RANDOM  ARRAYS 

As  discussed  previously,  the  problems  associated  with  the  implemen- 
tation of  high-resolution  image  systems  (whether  full  arrays,  uniformly 
thinned  arrays,  or  algorithmically  designed  thinned  arrays)  may  be 
alleviated  by  using  thinned  random  arrays.  Random  arrays  offer  distinct 
advantages  over  conventional  arrays  not  only  in  their  ability  to 
accomplish  the  objectives  of  the  imaging  system  with  many  fewer  elements, 
but,  perhaps  equally  important,  also  in  their  ability  to  predict  the 
a priori  probability  of  a successful  design.  The  mathematical  details 
of  the  properties  of  random  arrays  have  been  investigated  by  Lo  as  they 
apply  to  antenna  theory.  Some  essential  results  of  this  work  are  repro- 
duced in  Appendix  A.  Although  this  theory  has  been  well  developed  for 
antenna  applications,  its  usage  in  imaging  systems  has  not  been  explored 
(as  evidenced  by  a search  of  the  published  literature).  The  many  attrac- 
tive advantages  of  this  class  of  thinned  arrays  are  discussed  below  as 
they  apply  to  high-resolution  imaging  systems. 
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Random  arrays  can  achieve  high  resolution,  low  sidelobes,  and  high 
directive  gain  with  only  a small  fraction  of  the  number  of  array  ele- 
ments normally  required  by  uniform  arrays  to  accomplish  the  same  objec- 
tives. Furthermore,  the  fraction  of  reduction  possible  (i.e.,  thinning 
factor)  becomes  more  dramatic  as  the  number  of  elements  normally  required 
in  the  uniform  array  increases.  For  this  class  of  arrays,  the  sidelobe 
levels  are  closely  related  to  the  number  of  elements  used  in  the  array 
and,  to  a much  lesser  degree,  to  the  aperture  dimension.  Steinberg  has 
concluded  that  the  peak  sidelobe  level  is  linearly  dependent  on  N,  the 
number  of  array  elements,  and  only  logarithmically  dependent  on  the 
aperture  dimension.  Therefore,  for  a given  number  of  elements,  higher 
and  higher  resolutions  can  be  obtained  by  spreading  these  elements  over 
a larger  and  larger  aperture;  the  sidelobe  levels  would  remain 
substantially  the  same  while  the  mainlobe  is  narrowed  proportionately. 

On  the  other  hand,  if  the  array  aperture  size  is  fixed,  the  required 
specifications  of  the  array  can  be  attained  with  many  fewer  elements  by 
using  random  arrays.  As  few  as  15  elements  will  give  good  agreement  to 
theory  under  c-  cain  conditions. 

Equally  important,  the  question  of  the  successful  realization  of 
a design  to  prescribed  specifications  must  be  considered.  Since  no 
general  theory  exists  for  the  algorithmic  design  of  aperiodic  arrays  to 
meet  prescribed  performance  levels,  the  results  must  necessarily  be 
obtained  by  a trial-and-error  procedure,  even  with  the  aid  of  high-speed 
computers.  Contrast  this  with  random  arrays,  which  give  an  a priori 
probability  of  success  in  meeting  design  specifications.  The  probability 
of  success  is  determined  before  any  detailed  computations  are  carried 
out. 

Although  the  properties  of  random  arrays  are  approached  from  a 
probabilistic  point  of  view,  the  performance  of  the  array  is  not 
probabilistic.  The  performance  of  the  array  is  completely  deterministic 
once  the  disposition  of  the  array  elements  has  been  decided.  Thus,  in 
this  approach,  the  design  of  an  array  is  reduced  to  playing  a game  of 
chance,  for  which  the  odds  in  favor  of  success  could  be  designed  to  be 
overwhelmingly  in  our  favor.  For  one-dimensional  arrays,  experimental 
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verification  of  these  characteristics  have  been  made  cn  antenna  patterns 
in  the  microwave  region  by  Lo  and  Simcoe.  The  agreement  between  experi- 
ment and  theory  was  remarkable.  Similarly,  Steinberg  has  compared  the 
peak  sidelobe  levels  between  random  arrays  and  algorithmically  designed 
aperiodic  arrays.  Specifically,  he  compared  70  algorithmic  arrays  with 
170  random  arrays  for  the  distribution  of  their  peak  sidelobes.  He 
found  that  both  distributions  were  nearly  log  normal  with  the  same 
average  and  median  values.  However,  a crucial  difference  existed  in 
their  standard  deviations.  The  standard  deviation  of  the  random  array's 
distribution  was  found  to  be  approximately  half  that  of  the  algorithmic 
array's.  Thus,  "the  compactness  of  the  raidom  distribution  almost 
guarantees  against  selection  of  a random  ar.*ay  with  catastrophically 
large  peak  sidelobes." 

Other  significant  advantages  exist  for  random  arrays.  Because  many 
fewer  elements  are  required  in  one  mode  of  application  of  random  arrays, 
the  rate  of  transmission  of  information  from  such  arrays  would  be  less 
than  that  required  for  full  or  uniform  arrays,  thereby  reducing  the 
bandwidth  of  the  transmission.  The  tolerances  on  the  exact  disposi  ion 
of  the  array  elements  need  not  be  very  stringent  for  random  arrays  to 
achieve  the  design  specifications.  Recognizing  the  fact  that  the  sidelobes 
of  random  arrays  are  disposed  at  random  locations  relative  to  the  mainlobe, 
there  exists  a very  simple  scheme  of  image  enhancement  to  suppress  the 
spurious  background  levels  caused  by  the  sidelobes.  By  averaging  the 
reconstructed  images  obtained  from  different  random  array  configurations, 
the  background  levels  can  be  suppressed  as  U/ft,  where  M is  the  number  of 
images  used  for  the  averaging  process.  The  different  array  configurations 
can  be  achieved  through  the  use  of  reflector  elements  placed  at  various 
locations,  and  the  signals  from  these  reflectors  can  be  transmitted  to  the 
receiver  elements  through  a switching  sequence.  Finally,  by  the  very 
nature  of  the  random  array  concept,  the  performance  of  these  arrays  is 
not  very  susceptible  to  individual  array  element  failures  or  destruction. 
Subset^  of  larger  random-array  ,'onf  igurations  will  perform  adequately 
without  catastrophic  effects.  This  graceful  degradation  is  shown  in  the 
computer  simulation  results  in  Section  2.D. 


The  incoherent  OTF  uf  the  thin  random  array,  which  can  be  calculated 
by  its  autocorrelation,  will  have  very  few  redundant  vector  spacings 
(see  Figures  5,  7,  and  8).  There  are  N(N-l)/2  vector  spacings  produced 
by  N array  elements.  For  .ncoherent  imaging,  these  spacings  in  autocor- 
relation space  are  the  spatial  frequency  sampling  points.  The  OTF 
measures  the  existence  and  number  or  redundancy  of  these  samples  as  a 
function  of  spatial  frequency.  If  it  were  possible  to  sample  fully  and 
have  no  redundancy  at  all,  then  it  would  require  approximately  only  ^2N 
array  elements  to  have  N desired  unique  samples.  The  full  array  of  N 
elements  can  distinguish  N unique  resolution  elements  in  its  FOV.  This 
may  be  seen  by  noting  that  the  FOV  of  the  array  is  determined  by  the 
diffraction  angle  of  the  array  element  and  is  inversely  proportional  to 
the  element  dimension  while  resolution  is  proportional  to  array  dimen- 
sion. For  example,  for  a 4096-element  full  array,  only  /2N  'V*  91  elements 
or  /2/N  'v  2%  of  the  full  array  would  be  required  for  complete  sampling 
were  it  not  for  occasional  redundancies  and  missing  spatial  frequencies. 
More  elements  are  needed  to  make  up  the  deficiency.  In  one  dimension, 
there  are  only  a very  few  nonredundant  arrays  possible.  In  two  dimen- 
sions, it  is  impossible,  in  principle,  to  have  a nonredundant  array  that 
fills  the  autocorrelation  space  of  the  array.  In  our  studies  for  mono- 
chromatic imaging,  we  have  found  it  necessary  to  employ  random  array 
elements  only  three  times  more  numerous  than  this  theoretically  ■ 
unattainable  nonredundant  design  number.  For  polychromatic  imaging  (as 
explained  below),  only  1.5  times  more  was  needed.  Thus,  thinned  random 
arrays  can  automatically  perform  low  redundancy  sampling. 

The  above  discussion  has  concerned  monochromatic  or  quasi- 
monochromatic  radiation.  For  polychromatic  incoherent  imaging,  even 
fewer  elements  are  needed  in  the  thinner  array  to  get  a good  sampling  of 
the  spatial  frequency  plane.  Each  array  element  pair  defines  or  measures 
a different  spatial  frequency  for  each  different  wavelength, so  many  new 
spatial  frequencies  are  sampled  as  the  bandwidth  is  increased.  If  the 
object  does  not  have  a pathologic  spatial  frequency-"color"  relationship, 
polychromatic  imaging  with  a thin  random  array  wili  be  superior  to 
monochromatic.  Our  computer  simulations  bear  this  out  strikingly. 
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Another  class  of  arrays  may  be  specially  constructed  to  be  minimally 
redundant  or  even  nonredundant  and  compact  up  to  a maximum  spatial  frequency 
that  is,  however,  much  less  than  that  given  by  the  maximum  dimension  of  the 
array  (see  Figure  10).  This  class  of  arrays  shares  some  properties  with 
random  arrays,  but,  since  they  are  generated  algorithmically,  it  is  not 
proper  to  call  them  random.  These  two  classes  of  arrays  can  be  relatively 
immune,  in  a special  sense,  to  phase  aberrations. 

To  see  how  this  immunity  may  come  about,  consider  a pair  of  aperture 
elements  at  a fixed  spacing  that  defines  a certain  spatial  frequency,  f, 
as  the  elements  of  an  interferometer  producing  fringes  of  a certain 
"visibility."  The  MTF  of  the  aperture  is  a synthesis  of  the  visibilities 
of  all  such  pairs  for  all  frequencies.  Adding  an  additional  redundant 
pair  at  a frequency  f,  in  the  absence  of  phase  distortion,  always 
enhances  the  MTF  at  that  frequency.  But  consider  an  additional  pair  at 
f translated  somewhere  in  the  aperture  where  it  suff  rs  a relative 
space  variant  phase  distortion  due  to  the  atmosphere  or  distortion  in 
figure.  It  still,  as  a pair,  produces  fringes  of  the  same  visibility 
although  the  fringes  are  shifted.  That  is,  they  are  shifted  in  phase. 
Synthesizing  this  pair  with  the  original  pair  now  produces  a lesser 
resultant  in  the  MTF  than  it  would  in  the  absence  of  distortion  and 
produces  possibly  even  a reduction  over  that  produced  by  the  single  pair. 

On  the  average,  assuming  random  phases,  the  growth  of  the  MTF,  before 
normalization,  for  N contributing  pairs  of  elements  would  be  as  in  random 
walk  (i.e.,  proportional  to  ✓'n)  . But,  since  the  normalization  at  zero 
frequency  reduces  the  entire  transfer  function  by  1/N,  the  net  result  is 
a l/v'ft  reduction  in  the  MTF  at  that  spatial  frequency.  For  nonrandom 
phase  disturbances,  the  MTF  may  become  very  3mall,  or,  in  some  cases, 
negative.  Clearly,  nonredundant  or  minimally  redundant  arrays  would  not 
suffer  significantly  from  these  problems.  Radio  astronomers  have  been 
interested  in  this  possibility.  It  remains  to  be  decided  by  experiment 
whether,  in  the  presence  of  space  variant  phase  distortion,  it  is  more 
desirable,  with  respect  to  given  image  quality  criteria  for  given 
classes  of  images,  to  have  the  arbitrary  sampling  of  space  varying  phase 
in  the  thin  array  rather  than  to  have  tue  average  but  distorted  phase  in 
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the  full  array.  One  such  comparison  experiment  is  described  in 
Section  2.E.2. 

D.  SIMULATION  PROCEDURE 

Computer  simulation  studies  of  random  array  imaging  were  carried 
out  for  the  incoherent,  coherent,  polychromatic  incoherent,  and  parrialy 
coherent  cases.  A modular  software  package  was  designed  for  use  on  the 
simulation  studies.  This  approach  was  taken  to  provide  a flexible  single 
executive  control  program  that  would  simulate  these  cases  while  allowing 
optional  system  parameters  to  be  introduced  as  needed.  These  optional 
system  parameters  include  Gaussian  or  Poisson  noise,  atmospheric  turbulence, 
aperture  configuration,  and  percent  thinning  factor  in  the  thinned  random 
arrays.  Additionally,  the  output  format  of  the  reconstructed  image  and 
pertinent  intermediate  data  files  can  be  chosen  to  be  in  floating  point 
or  packed  integer  arrays ^ the  latter  with  or  without  a header.  A block 
diagram  of  the  simulation  software  package  is  shown  in  Figure  3. 

The  simulations  are  based  on  the  Fourier  representations  of  the 
incoherent  and  coherent  imaging  process,  shown  by  Eqs.  4 and  8.  The 
basic  components  of  the  simulation  are  a Fourier  transform  of  the 
digitized  object  array  (i.e.,  the  input  image),  OTF  corresponding  to 
either  the  incoherent  or  coherent  case,  and  an  inverse  transform  of  the 
product  of  the  above  two  components  to  reconstruct  the  object  (i.e., 
the  reconstructed  image) . Additional  components  of  the  simulation  are 
incorporated,  as  needed,  to  determine  the  influence  of  system  parameters 
on  image  quality.  The  sequence  of  operations  that  takes  place  in  the 
simulations  is  described  below. 

For  the  coherent  imaging  case,  a random  phase  of  range  0 to  2tt  is 
first  multiplied  on  a pixel-to-pixel  basis  to  the  digitized  input  image. 
This  simulates  the  effect  of  speckle,  which  appears  in  coherently 
illuminated  objects.  The  transformation  of  the  input  image  to  the  far- 
field,  or  aperture,  plane  is  accomplished  by  taking  a fast  Fourier  trans- 
formation of  the  input  image.  At  this  point,  an  aperture  configuration 
is  generated  that  consists  of  a binary  pattern  of  regular  or  random  array 
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elements.  Examples  of  regular  and  random  array  patterns  are  shown  in 
Figure  4.  The  patterns  are  represented  in  binary  form,  with  bright  pixels 
corresponding  to  a value  of  1 and  dark  pixels  to  a value  of  0.  A value 
of  1 in  the  aperture  configuration  represents  a point  where  an  antenna 
element  in  the  array  exists,  and  a value  of  0 represents  the  absence  of 
an  antenna  element.  The  aperture  plane  spatial  frequencies  are  sampled 
only  at  those  locations  where  an  antenna  element  exists  (i.e.,  where  the 
antenna  configuration  is  represented  by  a 1 in  the  binary  array).  The 
product  of  the  Fourier  transform  of  the  input  image  with  the  aperture 
configuration,  taken  on  a pixel-by-pixel  basis,  then  produces  the  Fourier 
representation  of  the  reconstructed  image.  Finally,  an  inverse  Fourier 
transform  of  this  product  array  produces  the  reconstructed  image,  as 
imaged  with  the  corresponding  antenna  array  configuration. 

For  the  incoherent  imaging  case,  there  exist  two  major  differences 
in  the  simulation  procedure.  First,  there  is  no  complex  random  phase 
associated  with  the  input  image  because  an  incoherently  illuminated 
object  produces  only  an  intensity  distribution  across  the  object.  All 
phase  information  is  suppressed.  The  fast  Fourier  transform  is  taken 
with  respect  to  the  intensity  distribution  of  the  input  image  only. 

Second,  the  OTF  for  the  incoherer t.  case  is  no  longer  a binary  array 
representative  of  the  antenna  configuration  above.  For  incoherent 
imaging,  the  OTF  is  derived  by  taking  a Fourier  transform  of  the  aper- 
ture configuration,  then  a complex  absolute  value  squared  of  this 
Fourier  transform,  and  finally  a second  Fourier  transform  of  this  func- 
tion. In  coordinate  space,  this  is  equivalent  to  taking  an  auto- 
correlation of  the  aperture  configuration.  Thus,  the  spatial  frequencies 
in  the  aperture  plane  are  sampled  by  pairs  of  antenna  array  elements, 
with  the  sampled  spatial  frequencies  represented  by  vector  distances 
between  all  possible  pairs  of  elements  in  the  array.  Thus,  depending  on 
the  aperture  configuration,  redundant  sampling  is  likely  to  occur.  With 
the  exception  of  the  above  two  differences,  the  simulation  scheme  for 
coherent  and  incoherent  imaging  is  identical. 
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(a)  REGULAR  ARRAY 

64  x 64  APERTURE  SIZE 
4096  ELEMENTS 


(b)  RANDOM  ARRAY 

64  x 64  APERTURE  SIZE 
512  ELEMENTS  - 12.5% 


Figure  4.  Representative  aperture  configurations. 
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The  incoherent  MTF  of  a full  64  square  aperture  is  shown  in  Figure  5. 

Note  the  high  peak  (large  redundancy)  plotted  in  the  center  of  the  pattern 

at  zero  spatial  frequency  and  the  general  fall-off  until  the  aperture- 

limited  frequency  is  reached.  The  PSF  of  this  full  square  array  is  shown 

2 

in  Figure  6 and  exhibits  the  familiar  sine  behavior.  For  comparison, 

2 

the  MTF  of  a random  array  of  12.5%  of  the  64  elements  is  shown  in 

2 

Figure  7 (a  and  b).  Only  the  central  64  is  plotted  there.  Note  tl_at, 

except  for  the  central  peak  at  zero  frequency,  the  spatial  frequency 

plane  is  almost  fully  uniformly  sampled  with  the  low  redundancy. 

Although  not  shown  here,  the  sampling  continues  to  the  boundary.  The 

MTF  of  a 100-element  random  array  (2.44%)  is  shown  in  Figure  8 plotted 

2 

in  slightly  more  than  one  quadrant  (80  ) . Sampling  to  the  boundary  is 
evident  here.  The  PSF  of  the  12.5%  random  array  is  shown  in  Figure  9. 

It  has  a central  peak  width  almost  identical  to  the  central  peak  width  of 
the  full  array  and  many  additional  side  peaks  of  small  magnitude. 

Figure  10  shows  the  MTF  of  a small  25-element  algorithmically  designed 
array.  The  MTF  is  compact  and  nonredundant  up  to  a small  frac- 
tion of  the  frequency  domain  defined  by  the  full  aperture.  A more  com- 
plete sampling  is  impossible  in  principle.  For  minimal  redundancy  with 
almost  complete  sampling,  the  random  array  provides  an  effective  design 
cion.  An  illusion  makes  the  region  near  the  center  appear  to  be 
* \ os ’ . 

We  have  simulated  polychromatic  incoherent  imaging  by  the  process  of 
averaging  or  overlaying  several  monochromatic  images,  each  representing 
a difi  rent  wavelength.  In  the  simulation,  different  spectral  compo- 
nentrr  rere  modeled  by  different  overall  aperture  dimensions.  This  is 
legitimate  because,  in  the  far  field  of  diffraction,  wavelength  and 
aperture  size  scale  directly.  For  example,  if  a full  aperture  of  64  x 64 
elements  represents  4096  resolution  elements  at  a certain  wavelength,  then 
a smaller  size  32  x 32  aperture  correctly  represents  the  smaller  1024 
resolution  elements  that  would  obtain  from  the  original  64  x 64  aperture 
if  the  original  wavelength  were  made  twice  as  long.  The  same  physical 
aperture  at  two  different  wavelengths  can  be  moaeled  by  two  different 
computer  apertures  of  linear  dimensions  inverse  to  the  wavelengths. 
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This  is  because  the  physical  diffraction  angles  are  determined  by  the 
ratio  of  physical  aperture  to  wavelength.  If  everything  else  is  physically 
constant,  red  light,  for  example,  would  have  less  resolving  power  than  blue. 

The  sizes  of  the  individual  elements  in  the  various  apertures  we  used 
were  kept  constant  at  a single  element  and  positioned  at  the  same  relative 
location  while  the  size  of  the  total  aperture  was  varied.  This  correctly 
models  the  constancy  of  field  of  view  with  changing  wavelengths.  Each 
image  was  renormalized  to  8 bits  to  achieve  constant  antenna  efficiency  as 
a function  wavelength.  Uniform  spectral  illumination  was  assumed  in  this 
study. 

For  the  case  of  partially  coherent  imaging,  the  simulation  makes  ase 
of  the  specially  constructed  separable  coherence  functions  described  in 
Section  2. A. 3.  An  arbitrary  coherence  function  may  be  approximated  in 
this  fashion.  But  if  it  is  not  so  represented,  in  separable  form,  the 
general  four-dimensional  convolutions  or  Fourier  transforms  that  would 

2 

then  be  needed  would  far  exceed  the  computer  storage  available  for  128 
size  images.  The  computer  simulation  previously  used  for  coherent  and 
incoherent  imaging  was  modified  to  accept  our  separable  partially  coherent 
model.  A block  diagram  of  this  simulation  is  shown  In  Figure  11.  The 
procedure  is  to  input  an  image  of  128  x 128  pixels,  represented  by  their 
intensity  values.  The  square  root  of  each  intensity  valu>-  is  taken  and, 
when  speckle  modeling  is  desired,  converted  to  a complex  number  with  a 
random  phase  factor.  This  gives  us  a model  for  the  complex  disturbance 
image.  This  image  was  then  multiplied  by  each  term  in  the  separable 
representation  of  the  degree  of  coherence  function  y (Eq.  22),  and  the 
Fourier  transform  was  taken  of  each  resultant  term.  Appendix  C gives 
an  example  of  the  terms  of  such  a separable  representation  for  y.  Each 
term  is  then  multiplied  by  the  optical  transfer  function  for  the  aperture 
followed  by  the  inverse  transform.  The  resulting  image  terras  were  then 
squared  and  added  to  obtain  the  output  intensity  image.  Since  only  the 
image  intensities  are  of  interest,  no  cross  terras  are  computed,  and  thus 
neither  the  degree  of  coherence  nor  mutual  coherence  of  the  output  image 
was  computed.  By  this  means,  the  impossible  burden  of  extreme  storage 
needed  for  four-dimensional  computation  was  transformed  to  a tractable 
but  lengthy  and  expensive  one. 
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Figure  11.  Block  diagram  of  partially  coherent  simulation. 


The  aperture  configuration  used  in  the  simulations  is  generated 
straightforwardly.  A regular  array  is  generated  by  setting  the  antenna 
element  coordinates  of  a predetermined  regular  pattern  to  1 and  the 
remaining  array  positions  to  0.  Generally,  a square  aperture  array  con 
figuration  was  used.  For  the  random  array,  a random  number  generator  is 
used  to  produce  the  sequence  of  numbers  used  as  coordinates  of  the  random 
array  elements.  Each  pair  of  random  numbers  forms  the  (x,y)  coordinates 
of  an  array  element.  The  random  number  generator  is  normalized  so  that 
the  range  of  numbers  generated  falls  within  the  limits  of  the  aperture 
size  chosen. 

Additional  optional  system  parameters  are  available  to  include  in 
the  simulation  studies  to  determine  the  effect  of  these  parameters  on 
image  quality.  These  include  a Poisson  noise  generator,  a Gaussian 
noise  generator,  and  an  atmospheric  turbulence  generator.  The  noise 
generators  were  used  tc  determine  the  effect  of  S/N  on  image  quality  for 
both  regular  and  random  array  Imaging.  The  turbulence  generator  was  used 
to  determine  the  immunity,  if  any,  of  random  arrays  to  atmospheric 
turbulence.  Finally,  there  are  several  options  available  to  oueput  the 
reconstructed  image  on  disk  for  display.  Either  floating  point  format 
or  integer  format  is  available,  with  the  latter  form  normalized  to 
3-bits  over  che  dynamic  range  of  the  reconstructed  image  to  display  it 
on  gray-scale  displays.  Options  to  output  the  aperture  configutV^on 
and  OTF  on  disk  are  also  available. 

E.  RESULTS  OF  THE  SIMULATION  STUDIES  4 

Both  incoherent,  coherent,  and  partially  coherent  cases  were 
simulated  for  monochomatic  illumination.  Polychromatic  incoherent  imaging 
was  simulated  as  well.  In  these  studies,  the  input  image  (object)  used 
was  a transparency  digitized  to  8 bits  (24  dB)  of  dynamic  range  and  con- 
taining a high  background  intensity  level.  This  image  was  digitized  to  a 
128  x 128  pixel  array,  the  limitation  on  the  size  of  the  array  being 
imposed  primarily  by  the  core  space  available  in  the  computer.  The  size 
of  the  largest  full  aperture  was  chosen  to  be  64  x 64  elements  centered 
in  the  aperture  array.  With  this  full  array  size,  the  resolution  element 
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in  the  input  image  is  2 x 2 pixels  in  size.  A total  of  64  * 64  x 64  m 

128  x 128/2  x 2 4096  resolution  elements  are  therefore  produced  in  the 

full  array's  FOV.  From  128  to  2048  elements  were,  used  in  the  random  array 
simulation  studies,  corresponding  to  3%  and  50%,  respectively,  of  the 
number  of  full  array  elements.  All  random  arrays  generated  in  the  studies 
were  confined  within  the  64  x 64  aperture  boundary  of  the  fail  array. 

Figures  12  through  16  show  the  results  of  incoherent  imaging  with 
thinned  random  arrays.  Included  in  each  figure  are  the  object  and  the 
full  aperture  diffraction-limited  image.  These  are  shown  for  comparison 
with  the  random  array  images.  Figures  17  through  18  show  the  results  of 
coherent  Imaging  with  thinned  random  arrays;  the  corresponding  full  array 
image  is  ihcluded  for  comparison.  Because  of  the  onoet  of  speckle  effects 
in  coherent  imaging,  an  additional  parameter,  the  number  of  overlays, 
is  introduced.  The  effects  of  speckle  can  be  reduced  by  speckle  averaging, 
in  which  reconstructed  coherent  images  with  different  speckle  patterns 
are  added  to  "smooth  out"  the  individual  speckle  patterns.  The  coherent 
images  are  shown  as  a function  of  the  number  of  overlays  used  for  speckle 
averaging. 

Figure  19  compares  full  and  random  arrays  for  incoherent  imaging  in 
the  presence  of  turbulence.  Figures  20  through  25  3how  the  results  of 
polychromatic incoherent  imaging  with  thinned  random  arrays  and  with  the 
corresponding  full  array  images  for  comparison.  Partially  coherent 
imaging  was  successfully  simulated  by  modeling  the  coherence  in  terms  of 
totally  separable  functions  to  avoid  four-dimensional  transforms,  which 
would  have  exceeded  the  storage  capabilities  of  our  computer.  Never- 
theless, tne  large  computational  burden  precluded  more  than  a cursory 
examination  of  the  problem.  A single  partially  coherent  image  required 
1 hr  of  CPU  on  a PDP  KL-10  with  a Tenex  operating  system.  Partially 
coherent,  fully  coherent,  and  incoherent  images  are  shown  in  Figure  26. 

1 . Incoherent  Imaging 

The  image  quality  achievable  with  random  arrays  thinned  to 
only  12.5%  of  the  number  of  full  aperture  elements  is  shown  in  Figure  12. 
Four  different  configurations  of  random  arrays  were  arbitrarily  chosen, 
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Figure  19.  Incoherent  imaging  in  presence  of  atmospheric  turbulence 
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ORIGINAL  OBJECT  MONOCHROMATIC  IMAGE  FROM 

RANDOM  ARRAY  NO.  2 

26%  OF  FULL  APERTURE  IMAGE 


Figure  20.  Incoherent  monochromatic  imaging. 
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MONOCHROMATIC  DIFFRACTION  LIMITED 
FULL  APERTURE  IMAGE 


RANDOM  ARRAY  NO.  1 

26%  OF  FULL  APERTURE  IMAGE 


FULL  APERTURE  IMAGE 


RANDOM  ARRAY  NO.  2 
26%  OF  FULL  APERTURE 


Figure  21.  Incoherent  polychromatic  imaging, bandwidth 
span  factor  1.14. 
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MONOCHROMATIC  DIFFRACTION  LIMITED 
FULL  APERTURE  IMAGE 


RANDOM  ARRAY  NO.  1 

25%  OF  FULL  APERTURE  IMAGE 
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FULL  APERTURE  IMAGE 


RANDOM  ARRAY  NO.  2 

25%  OF  FULL  APERTURE  IMAGE 


Figure  23. 


Incoherent  polychromatic  imaging,  bandwidth 
span  factor  2.00. 
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c)  MONOCHROMATIC  IMAGE  RANDOM  d)  POLYCHROMATIC  IMAGE  RANDOM 

ARRAY  12.5  % OF  FULL  APERTURE  ARRAY  12.5  % OF  FULL  ARRAY 

BANDWIDTH  FACTOR  TWO 


Figure  24.  Polychromatic  incoherent  imaging,  100%  and  12.5%  of  full 
array,  bandwidth  span  factor  2.00. 
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MONOCHROMATIC  IMAGE  RANDOM 
ARRAY  6.25  % OF  FULL  APERTURE 
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b)  POLYCHROMATIC  IMAGE  RANDOM 
ARRAY  6.25%  OF  FULL  APERTURE, 
BANDWIDTH  FACTOR  TWO 


c)  MONOCHROMATIC  IMAGE  RANDOM  d)  POLYCHROMATIC  IMAGE  RANDOM 

ARRAY  3.125  % OF  FULL  APERTURE  ARRAY  3.125  % OF  FULL  APERTURE 

BANDWIDTH  FACTOR  TWO 


Figure  25.  Polychromatic,  incoherent  imaging,  6.25%  and  3.125% 
of  full  array,  bandwidth  span  factor  2.00. 


a)  COHERENT  64  FULL  ARRAY 


Figure  26.  Partially  coherent  imaging 


c)  PARTIALL  COHERENT  642 
FULL  ARRAY 


Figure  62.  (Continued) 


each  containing  512  elements.  The  full  ax ray  is  composed  of  4096  elements. 
The  four  aperture  configurations  were  generated  by  using  the  first  four 
sequences  of  random  numbers  produced  by  the  random  number  generating  rou- 
tine. High  signal-to-noise  ratios  were  used  in  this  simulation  to  isolate 
the  thinning  factor  effects  from  the  noise  effects  in  degrading  image 
quality.  The  image  quality  achievable  with  random  arrays  thinned  to  12.5% 
of  the  full  array  compares  favorably  with  the  quality  achievable  with  the 
full  diffraction-limited  array. 

One  advantage  of  random  arrays  is  that  the  image  quality  does  not 
suffer  catastrophic  degradation  when  a large  fraction  of  the  array  ele- 
ments fails  (see  Section  1).  The  graceful  degradation  of  random  array 
images  when  large  fractions  of  array  elements  are  removed  is  shown  in 
Figure  13.  Here,  the  random  array  elements  were  removed  randomly,  with 
50%  of  the  elements  removed  each  time.  Starting  at  1024  elements  in  a 
random  array  thinned  to  25%  of  the  full  array  number.  Figure  13(b,c,e,f) 
represents  random  arrays  thinned  to  25%,  12.5%,  6.25%,  and  3.125%  of  the 
full  array.  In  each  case,  when  random  arrays  with  fewer  elements  are 
created,  these  arrays  are  subsets  of  the  previous  larger  arrays.  An 
examination  of  the  random  array  images  in  Figure  13  shows  that  there  is 
no  drastic  change  in  image  quality  until  thinning  has  been  reduced  to 
below  the  5%  level  (Figure  13(f)).  High  S/N  conditions  also  prevailed 
in  this  sequence  of  simulations. 

Another  advantage  of  thinned  random  arrays  (one  mentioned  previously) 
is  that  system  resolution  can  be  increased  dramatically  with  the  same 
number  of  array  elements  in  use.  All  that  is  necessary  is  to  spread  the 
fixed  number  of  elements  over  a larger  aperture  area.  This  is  vividly 
depicted  in  Figure  14.  A 400-element  full  array  was  used  as  a reference 
for  comparison.  The  resolution  achievable  with  a full  400-element  array, 
with  aperture  size  of  20  x 20  pixels,  is  shown  in  Figure  14(a).  Merely 
redistributing  the  400  elements  onto  larger  size  aperture  limits,  as  was 
done  in  Figure  14(b,c),  significantly  improved  resolution  without  notice- 
ably degrading  image  quality.  Here,  a threefold  increase  in  resolution 
was  achieved  (see  Figure  14 (a  and  c)).  For  a comparison  of  image  quality, 
full  array  diffraction-limited  images  corresponding  to  the  increased 
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aperture  sizes  are  shown  in  Figure  14(e,f).  These  two  full  arrays  were 
comprised  of  1936  and  4096  elements,  respectively. 

Figures  15  and  16  show  the  effect  of  low  S/N  ratio  conditions  on  image 
quality.  A Gaussian  noise  spectrum  was  added  to  each  array  element  to 
simulate  the  effect  of  receiver  noise.  The  variance  of  the  Gaussian 
spectrum  was  adjusted  relative  to  the  average  intensity  of  the  Fourier 
components  to  produce  the  S/N  ratios  shown.  Figures  12 (c,  d)  and 
13 (c,  d)  represent  the  image  quality  obtained  from  random  arrays  with 
512  elements  (12. 5%  of  the  full  array)  for  the  S/N  ratios  shown.  The 
corresponding  full  array  im:  or  identical  S/N  ratios  are  shown  below 

each  random  array  image.  For  \'/y.  ratios  above  3,  the  random  array  image 

quality  differs  little  £?.-.«■ m the  lull  array  Images. 

2.  Coherent  Imaging 

Results  of  the  computer  simulation  studies  of  coherent  imaging  with 
random  and  regular  arrays  are  shown  in  Figures  17  and  18.  In  this 
sequence  of  simulations,  the  number  of  overlays  is  varied  to  determine 
the  effect  of  speckle  averaging  on  image  quality.  In  each  figure,  the 
random  array  Images  are  shown  on  the  right  side,  with  the  corresponding 
full  array  images  shown  on  the  left.  An  aperture  size  of  64  x 64  pixels 

is  used,  and  the  random  array  was  thinned  to  50%  of  the  number  of  full 

array  elements.  A high  S/N  ratio  prevailed  in  these  simulations. 

For  the  coherent  Imaging  case,  the  OTF  is  the  aperture  configuration 
itself.  Therefore*  when  an  element  of  the  array  is  removed,  as  in  a 
thinned  array,  the  spatial  frequency  corresponding  to  that  array  element 
Is  not  sampled  at  all. 

Speckle  effects  severly  degrade  image  quality.  For  the  random 
array  used,  speckle  effects  still  are  prominent  even  after  100  overlays, 
while  for  the  full  array,  speckles  can  be  significantly  removed  at  100 
overlays.  For  lower  numbers  of  overlays  (e.g.,  less  than  10  in  this 
case),  the  speckles  are  quite  severe  in  both  situations. 


3.  Turbulence  Effects 


The  effect  of  atmospheric  turbulence  on  image  quality  obtainable 
from  random  and  regular  array  incoherent  imaging  is  shown  in  Figure  19 
for  which  a random  array  thinned  to  12.5%  of  the  full  array  elements  was 
used.  Turbulence  levels  ranging  from  moderate  to  severe  were  used*  The 
random  array  images  are  again  shown  on  the  left  column  and  the  full  array 
images  on  the  right.  The  image  quality  obtained  with  the  12.5%  random 
arrays  is  degraded  more  severely  than  the  corresponding  full  array  images. 
For  moderate  turbulence  levels  (Cat  0.5),  the  difference  in  image  quality 
is  slight  between  the  random  and  full  arrays  for  higher  levels  of  turbu- 
lence, but  since  images  are  severely  degraded  anyway,  comparisons  between 
them  are  less  meaningful.  There  seems  to  be  no  evidence  from  this  experi- 
ment to  support  the  idea  that  the  phase  aberration  immunity  for  thinned 
random  arrays,  discussed  in  Section  2.C,  is  of  any  value  for  images  of 
this  class. 

4.  Polychromatic  Incoherent  Imaging 

The  results  of  our  computer  simulation  studies  of  polychromatic 
incoherent  imaging  with  random  arrays  are  summarized  below.  Typical 
results  are  shown  in  Figures  20  through  25.  In  these  studies,  the  object 
used  was  a digitized  photograph  with  8 bits  (24  dB)  of  dynamic  range;  it 
contained  a high  background  intensity  level.  The  image  was  digitized  to 
a 128  x 128  pixel  array.  The  size  of  the  full  aperture  for  the  short 
wavelength  limit  of  the  polychromatic  bandwidth  was  chosen  to  be  64  x 64 
elements,  thereby  producing  4096  resolution  elements  in  the  system  FOV  at 
this  wavelength.  In  every  polychromatic  composite  average,  the  percentage 
of  the  number  of  elements  of  each  monochromatic  component  is  kept  fixed 
rather  than  the  number  of  elements.  The  number  of  elements  used  in  the 
random  array  simulation  studies  ranged  from  1024  to  32  elements,  cor- 
responding to,  in  separate  studies,  from  25%  to  3.125%  of  the  number  of 

2 2 

full  array  elements  for  apertures  varying  from  64  to  32  in  dimension. 

This  range  corresponded  to  a factor  of  two  in  wavelength.  Each  figure 

2 

(20  through  24)  includes  the  monochromatic  diffraction-limited  64  full 
aperture  image.  This  is  shown  as  a reference  for  comparison. 
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Figure  20  shows  two  images  found  with  two  different  random  arrays 

2 

each  thinned  to  25%  of  the  64  aperture,  giving  1024  elements.  These  are, 
of  course,  monochromatic  images  and  show  no  appreciable  degradation  due  to 
thinning. 

Figure  21  compares  the  polychromatic  images  found  by  averaging  full 
9 2 2 

arrays  of  64  and  56  elements  with  those  from  two  different  25%  thinned 
arrays  of  the  same  dimensions.  The  effect  of  combining  the  two  simulated 
wavelengths  is  very  small,  and  the  effects  of  thinning  are  entirely  neg- 
ligible. In  Figure  22  similar  results  are  seen  for  arrays  of  dimension 

64^,  56^,  and  48^  each  thinned  to  25%  of  64^,  56^,  and  48^,  respectively. 

2 2 2 2 2 
In  Figure  23  for  arrays  of  dimension  64  , 56  , 48  , 40  , and  32  , 

which  represent  a span  of  a factor  of  two  in  wavelength  or  bandwidth, 

the  effects  of  averaging  the  more  blurry  long  wavelength  image  with  the 

sharper  short  wavelength  image  are  more  noticeable.  But  the  effects  of 

thinning  the  full  array  are  still  negligible. 

2 

Figure  24  again  shows  a monochromatic  diffraction-limited  64  full 

aperture  image  as  a new  reference  for  comparing  this  and  the  next  figure. 

Figure  25.  In  Figure  24(a),  a more  continuous  polychromatic  image  is 

shown;  it  was  generated  from  10  different  monochromatic  full  array  over- 

2 

lays  spanning  a bandwidth  factor  of  two  in  the  array  sizes  of  from  64  to 
2 

32  . A slight  loss  of  resolution  may  be  noted  in  comparison  to  the  short 

wavelength  monochromatic  image  in  Figure  24(a).  In  Figure  24(c),  the 

monochromatic  image  was  produced  with  a random  12.5%  of  the  short- 
2 

wavelength  (64  ) aperture.  Some  small  degradation  in  comparison  with 

the  full-array  image  can  be  noted.  The  polychromatic  image.  Figure  24(d), 

although  showing  less  resolution  than  Figure  24(c),  does  not  show,  for 

example,  the  ghosting  evident  in  the  background. 

A similar  comparison  is  made  in  Figure  25 (a  and  b)  for  thinning  to 

2 

6.25%.  The  degraded  short  wavelength  (64  ) monochromatic  image  of  Fig- 
ure 25(a)  is  compared  to  the  polychromatic  version  of  bandwidth  factor 
two  in  Figure  25(b).  The  severly  degraded  3.125%  short  wavelength  (64  ) 
monchromatic  image  of  Figure  25(c)  is  very  much  improved  in  the  polychromatic 
factor  two  bandwidth  version  in  Figure  25(d). 
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5.  Partially  Coherent  Imaging 


The  effects  of  partial  coherence  on  image  formation  were  studied 
by  employing  separable  monochromatic  partial  coherence  functions  of  the 
kind  illustrated  in  Figure  2.  In  one  set  of  experiments,  speckle 
formation  was  studied  as  a function  of  the  degree  of  coherence.  As  was 
expected,  the  variance  of  the  speckle  distribution  function  decreased  as 
the  coherence  was  reduced.  The  dramatic  difference  between  coherent  and 
incoherent  Imaging  may  be  seen  by  comparing  Figure  25(a)  with  Figure  26(b) 
respectively.  The  coherent  picture,  without  any  added  random  phase, 
displays  the  vivid  edge  "ringing"  expected  from  the  diffraction  of  sharp 
boundaries.  This  can  be  noted  in  the  border  of  the  picture  in  the  face 
of  the  subject,  and  along  the  right  side  of  the  subject's  hair.  The 
incoherent  image.  Figure  26(b),  shows  no  obvious  comparable  effects.  The 
partially  coherent  image.  Figure  26(c),  using  the  121-term  separable 
coherent  function  plotted  in  Figure  2,  shows  an  intermediate  degree  of 
diffraction  effects.  These  effects  properly  are  somewhat  exaggerated 
when  discrete  as  opposed  to  continuous  arrays  are  employed  (as  noted  in 
Appendix  B) . 
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SECTION  3 


SUMMARY  OF  RESULTS 


The  results  from  the  computer  simulation  studies  of  random  array 
imaging  suggest  that  random  arrays  are  primarily  useful  for  incoherent 
imaging  applications.  The  advantages  of  random  array  imaging  for  inco- 
herent Imaging  are 

Large  degree  of  thinning  possible.  In  this  study,  we  have  shown 
that  random  arrays  thinned  to  less  than  6X  of  the  number  of 
full  array  elements  (4096)  produce  Incoherent  image  quality 
comparable  to  the  full  array  diffraction-limited  image  quality. 
The  degree  of  thinning  allowable  is,  in  principle,  inversely 
proportional  to  the  square  root  of  the  number  of  elements  in 
the  full  array,  which  number  is  equivalant  to  the  number  of 
elements  resolved  by  the  full  array  in  its  field  of  view. 

Polychromatic  radiation.  Thinning  to  3%  was  shown  to  be 
similarly  effective  in  the  case  of  broadband:  with  polychromatic 
imaging. 

Easy  to  design.  Many  random  array  configurations  produce  image 
qualities  comparable  to  full  array  diffraction-limited  quality. 
The  probability  of  success  can  be  estimated  for  choosing  an 
array. 

Graceful  degradation.  The  image  quality  obtainable  from 
random  arrays  does  not  suffer  catastrophic  degradation  when 
large  fractions  of  array  elements  are  removed.  The  image 
quality  degrades  very  gracefully. 

Increased  resolution  with  fixed  number  of  array  elements. 

The  resolving  capability  of  random  arrays  can  be  increased 
significantly  using  the  same  number  of  array  elements.  This 
la  accomplished  by  merely  redistributing  the  array  elements 
over  a larger  aperture  size. 


Partially  coherent  imaging  proved  extraordinarily  expensive  to 
simulate,  but  a separable  version  of  the  coherence  function  was  modeled 
that  allowed  a cursory  examination  of  the  problem.  As  expected,  the 
behavior  of  partially  coherent  images  was  qualitatively  between  the 
coherent  and  incoherent  extremes. 
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For  incoherent  Imaging  in  turbulent  conditions  in  the  atmosphere, 
thinned  random  arrays  suifer  only  slightly  in  image  quality  when  compared 
to  full  array  images  for  a large  range  of  turbulence  levels.  For  lower 
turbulence  levels,  the  difference  in  image  quality  is  slight,  although 
the  random  array  images  were  more  degraded.  For  higher  turbulence 
levels,  both  arrays  suffered  severe  degradations,  and  the  comparison  of 
image  quality  under  those  circumstances  is  less  meaningful. 

For  coherent  imaging  applications,  random  arrays  do  not  offer  sig- 
nificant advantages  over  full  arrays.  This  is  because  speckle  effects  can- 
not be  reduced  enough  by  speckle  averaging  (or  overlays) , and  substantial 
degradation  of  imagv  quality  still  remains  after  numerous  overlays. 
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APPENDIX  A 


RANDOM  ARRAY  DESIGN  THEORY 

There  is  a body  of  microwave  antenna  lite' ■>  're  dating  from  the 
1960s  describing  the  theoretical  expectations  large  antenna  arrays 
with  randomly  spaced  elements.  The  tradeoff o between  resolution,  gain, 
and  sidelobe  level  have  been  analyzed  in  probabilistic  terms.  These 
ideas  have  been  experimentally  tested,  and  the  results  were  favorable. 
Comparisons  were  also  made  between  relative  efficacy  of  antennas 
designed  algorithmically  and  those  designed  by  random  placement  of  the 
elements;  the  latter  showed  remarkable  advantages.  This  appendix 
briefly  summarizes  some  important  points  from  the  work  of  Lo. 

The  sidelobe  level  is  related  to  the  number  of  elements  and  only 
slightly  to  the  aperture  dimension.  High  resolution  can  be  obtained 
with  few  elements,  and,  for  a given  number  of  elements,  high  resolution 
can  be  obtained  by  spreading  these  elements  without  a substantial  change 
in  the  sidelobe  level  or  the  directive  gain. 

The  probability  of  a sidelobe  level  below  r is  [l  * exp (-Nr2)]  ^a, 
where  N is  the  total  number  of  elements,  a is  the  aperture  dimension 
measured  in  wavelengths,  and  r is  the  sidelobe  level  with  the  mainbeam 
normalized.  Figure  A-l  is  a universal  plot  of: this  relationship  plot- 
ting probability  versus  rv^  for  a variety  of  aperture  dimensions  indexed 
by  q - log jq  a.  For  planar  arrays,  q is  equal  to  the  sum  of  q,  and 
This  family  tells  the  probability  of  finding  ("designing")  a random 
array  with  a specified  maximum  sidelobe  level  for  a given  aperture  and 
number  of  elements.  Figure  A-2  demonstrates  the  nearly  perfect  agree- 
ment between  the  probabilistic  theory  and  a computer  simulation  with 
randomly  positioned  antenna  arrays. 


Y.T.  Lo  and  R.J.  Simcoe,  IEEE  Trans.  Ant.  and  Prop.,  Vol.  Ap-15,  No.  2, 
March  1967,  pp.  231-235. 


CUMULATIVE  PROBABILITY 


Figure  A-l.  Theoretical  distribution  curves  for  r^N,  where  N * total 
number  of  elements,  and  r « sidelobe  level.  In  the  case 
of  a linear  array,  q - log^Q  (array  length  in  wavelengths); 
in  the  case  of  a planar  array,  q - log^Q  (aperture  area 
in  wavelengths  squared). 
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SIDELOBE  LEVEL,  dB 

Figure  A-2. 

Comparison  of  the  theoret'cal  and  experimental 
distributions  for  the  sidelobe  levels  of  linear 
arrays.  The  experimental  distribution  is 
obtained  from  180  pseudo-random  linear  arrays 
generated  by  two  planar  arrays, 
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APPENDIX  B 


DIFFRACTION  EFFECTS  IN  DISCRETE  COHERENT  IMAGES 


The  general  effects  of  diffraction  in  coherent  images  are  well 

known  but  there  are  novel  striking  effects  when  the  arrays  are  dLscrete 

rather  than  corttinuous.  Especially  strong  symmetric  effects  can  be 

seen  when  the  imaging  array  is  full  and  the  number  of  elements  in  the 

array  and  in  tHfe  object  is  commensurate.  In  Figure  B-l,  the  coherent 

1 

images  of  a 64  bright  uniform  object  produced  by  four  different  full 

2 2 

arrays  are  shown.  In  Figure  B-l (a),  a 64  full  aperture  was  used,  and  16 

bright  crossed  lines  appear  in  the  image.  In  Figures  B-l(b),  B-l(c), 

2 2 2 

and  B-l(d),  7 t 8 , and  9 full  apertures  were  employed,  respectively . 

The  four  bright  corners  in  the  images  are  much  more  distinct  in  the  case 

2 2 
of  the  8 aperture, which  is,  of  course,  commensurate  with  64  . The 

2 

diffraction  effects  of  the  64  aperture  can  be  noted  in  the  coherent  and 
partially  coherent  images  in  Figure  26. 


a)  642 


b)  7l 


c)  8" 


d)  9^ 


Figure  B-l.  Diffraction  effects  in  coherent  imaging  with  discrete 
arrays.  Images  at  a 64z  uniform  object  with  full 
apertures  of  sizes:  a)  642,  b)  72,  c)  82,  d)  92. 
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APPENDIX  C 


SEPARABLE  REPRESENTATION  OF  TEE  DEGREE  OF  COHERENCE  FUNCTION 
Let  (X,Y)  be  the  pixel  location  of  the  vector  IT  in  Eq.  22. 


A - 126  x 2~9 
o 

*2  " 120  X 

2-9 

A.  ■ 10  x 2 
4 

■9 

A±  * 210  x 2-9 

A^  * 43  x 

2-9 

A5  - 2-9  in 

Eq. 

21. 

a * 2~6 

The 

M * 121  terms 

for 

each  f in 
n 

Eq.  22 

are: 

(1) 

A 

o 

(18) 

sin(2Y) 

• 

(2) 

cos(X) 

• 

/2A  A, 
o 1 

(19) 

sin(3Y) 

• 

^77 

(3) 

cos(2X) 

• 

/2A  A0 
o / 

(20) 

sin(4Y) 

• 

/2A  A, 
o 4 

(4) 

cos(3X) 

• 

/2  A A - 
o 3 

(21) 

sin(5Y) 

• 

/2A  A*. 

O D 

(5) 

cos(4X) 

* 

/2A  A. 
o 4 

(22) 

cos(X-Y) 

• 

(6) 

cos(5X) 

• 

/2A  A" 
o 3 

(23) 

cos(X-2Y) 

• 

/K?2 

(7) 

sin(X) 

• 

/2A  A- 
o 1 

(24) 

cos (X-3Y) 

• 

^A1A3 

(8) 

sin(2X) 

• 

Jik  A- 
o z 

(25) 

cos (X-4Y) 

• 

(9) 

sin(3X) 

• 

/2A  A» 
o 3 

(26) 

cos (X-5Y) 

• 

(10) 

sin(4X) 

• 

/2A  A, 
o 4 

(27) 

cos (2X-Y) 

• 

A2A1 

(11) 

sin(5X) 

• 

/2A  A. 
o 5 

(28) 

cos(2X-2Y) 

• 

-C 

(12) 

cos (Y) 

• 

✓2A  A. 
o 1 

(29) 

coa (2X-3Y) 

• 

(13) 

cos(2Y) 

• 

/2A  A_ 
o 2 

(30) 

cos(2X-4Y) 

• 

(14) 

cos(3Y) 

• 

/2A  A„ 
o 3 

(31) 

cos (2X-5Y) 

• 

*Vs 

(15) 

cos(4Y) 

• 

/2A“A. 
o 4 

(32) 

cos(3X-Y) 

• 

*Vl 

(16) 

cos(5Y) 

• 

/ITT 

o 5 

(33) 

ccs(3X-2Y) 

• 

v^A^A^ 

(17) 

sin(Y) 

• 

Jik'h, 
o 1 

(34) 

cos(3X-3Y) 

• 

4 

79 


(35) 

cos(3X-4Y) 

*^3A4 

(59) 

sin(3X-3Y) 

• 

/a3X3 

(36) 

cos(3X-5Y) 

* 

(60) 

sin(3X-4Y) 

• 

(37) 

cos(4X-Y) 

*^4A1 

(61) 

sin(3X-5Y) 

• 

>/a3a3 

(38) 

cos(4X-2Y) 

^A4A2~ 

(62) 

sin(4X-Y) 

• 

/S?T 

(39) 

cos(4X-3Y) 

* *71 

(63) 

sin(4X-2Y) 

• 

*^A4A2 

(40) 

cos(4X-4Y) 

*?4 

(64) 

sin(4X-3Y) 

• 

(41) 

cos(4X-5Y) 

*4A5 

(65) 

sin(4X-4Y) 

• 

A4A4 

(42) 

cos(5X-Y) 

* /A^A^ 

(66) 

sin(4X-5Y) 

• 

^7 5 

(43) 

cos(5X-2Y) 

* 

(67) 

sin(5X-Y) 

• 

*71 

(44) 

cos(5X-3Y) 

• A5A3 

(68) 

ain(5X-2Y> 

• 

/k^r2 

(45) 

cos(5X-4Y) 

* v'X5A4 

(69) 

sin(5X-3Y) 

• 

*%A3 

(46) 

cos(5X-5Y) 

‘ /*?5 

(70) 

sin(5X~4Y) 

• 

*71 

(47) 

sln(X-Y) 

(71) 

sin(5X-5Y) 

• 

*77 

(48) 

sln(X“2Y) 

* t/ A^2 

(72) 

cos (X+Y) 

• 

(49) 

sin(X-3Y) 

• / A^A^ 

(73) 

c.os(X+2Y) 

• 

/ A^2 

(50) 

sln(X-4Y) 

• ^ 

(74) 

cos (X+3Y) 

• 

✓Aia3 

(51) 

sin(X-5Y) 

(75) 

cos(X+4Y) 

• 

(52) 

sin(2X-Y) 

(76) 

cos (X+5Y) 

• 

■Vs 

(53) 

sln(2X-2Y) 

>/,A2A2 

(77) 

cos (2X+Y) 

• 

/A2A3 

(54) 

ain(2X-3Y) 

^A2A3 

(78) 

cos(2Y+2Y) 

• 

(53) 

sin(2X-4Y) 

* A2A4 

(79) 

cos(2X+3Y) 

* 

y/  A2  A^ 

(56) 

sin(2X-5Y) 

• Vk^k^ 

(80) 

cos(2X+4Y) 

• 

*^A2A4 

(57) 

sin(3X-Y) 

* */A3Al 

(81) 

cos (2X+5Y) 

• 

/ ^2A5 

(58) 

sln(3X-2Y) 

y/  A^2 

(82) 

cos(3X+Y) 

• 

/A3A1 

(83)  cos(3X+2Y) 

(84)  coa(3X+3Y) 

(85)  cos(3X+4Y) 

(86)  cos(3X+5Y) 

(87)  cos (4X+Y) 

(88)  cos(4X+2Y) 

(89)  cos(4X+3Y) 
(9U)  cos(4X+4Y) 

(91)  cos(4X-v-5Y) 

(92)  cos (5X+Y) 

(93)  cos(5X+2Y) 

(94)  cos(5X+3Y) 

(95)  cos(5X+4Y) 

(96)  coa(5X+5Y) 

(97)  sin(X+Y) 

(98)  sin(X+2Y) 
(99>  sin (X+3Y) 

(100)  sln(X+4Y) 

(101)  sln(X+5Y) 

(102)  sin(2X+Y) 


^A3A3 


/A4A1 

*^A4A2 

•v; 


/A^Aj 

1^A5A2 

*^5a3 

''Vs 

/a~a~ 

A1A1 


'VS 


(103)  sin(2X+2Y) 

(104)  sin(2X+3Y) 

(105)  ain(2X+4Y) 

(106)  aln(2X+5Y) 

(107)  ain(3X+Y) 

(108)  sin(3X+2Y) 

(109)  sin(3X+3Y) 

(110)  sin(3X+4Y) 

(111)  sln(3X+5Y) 

(112)  sin(4X+Y) 

(113)  aln(4X+2Y) 

(114)  sln(4X+3Y) 

(115)  sln(4X+4Y) 

(116)  sln(4X+5Y) 

(117)  sln(5X+Y) 

(118)  sin(5X+2Y) 

(119)  sin(5X+3Y) 

(120)  ain(5X+4Y) 

(121)  sln(5X+5Y) 


•*71 


^A3A3 


•*4*2 

^A4A4 


™np' 
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